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ABSTRACT 


It  is  shown  that  the  Ritz  projection  onto  spaces  of  piecewise  linear 

o  1 

finite  elements  is  bounded  in  the  Sobolev  space,  W  ,  for  2  <  p  <  •».  This 

p  —  — 

°  1  2 

implies  that  for  functions  in  W  O  W  the  error  in  approximation  behaves 

P  P 

1  O 

like  0(h)  in  W  ,  for  2  <  p  <  «®,  and  like  0(hz)  in  L  ,  for 
p  —  -  P 

2  <  p  <  “>.  In  all  these  cases  the  additional  logarithmic  factor  previously 


included  in  error  estimates  for  linear  finite  elements  does  not  occur. 
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SIGNIFICANCE  AND  EXPLANATION 
This  paper  concerns  error  estimates  for  methods  of  approximating  the 
solution  of  a  partial  differential  equation.  The  method  in  question  is  the 
so-called  "finite  element  method ^rtiich  was  developed  by  structural  engineers 
and  is  now  widely  used  in  all  branches  of  engineering.  The  paper  refines 
previously  derived  estimates  of  the  error  in  "maximum  norm,"  i.e.  the  maximum 
error  (as  opposed  to  an  average  error).  The  paper  settles  certain  technical 
questions  as  to  the  rate  of  convergence  of  the  finite  element  method  in  this 


norm. 
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SOME  OPTIMAL  ERROR  ESTIMATES  FOR  PIECEWISE 
LINEAR  FINITE  ELEMENT  APPROXIMATIONS 

Rolf  Rannacher*  and  Ridgway  Scott** 

1.  Introduction  and  Results. 

Let  ft  be  a  convex  polygonal  domain  in  1R  ,  and  let  it.  -  {K},  0  <  h  <  h  <  1. 

n  —  o 

be  finite  triangulations  of  n  such  that  the  usual  regularity  condition  is  satisfied: 
^  triangles  K  e  only  meet  in  entire  common  sides  or  in  vertices.  Each 

£?1.3g3le  K  e  \  contains  a  circle  of  radius  Cjh  and  is  contained  in  a  circle  c2h, 
where  the  constants  c^,c2  do  not  depend  on  K  or  h. 

Corresponding  to  ir^,  we  define  the  finite  dimensional  subspace  sh  C  W1  by 

Sh  "  ^vh  e  W»  :  vh  is  linear  °n  each  K  *  *h>  » 

and  the  Ritz  projection  R  :  -*•  S.  by 

h  2  h 


(l.D 


(VR^,^)  =  (Vu,V*h),  V  #h  t  Sh  . 


Here  Ld  and  ^ '  1  i  P  i  mew,  are  the 


Lebesgue  and  Sobolev  spaces  on 


,  respectively,  w  in  the  subspace 

m,p  p  *■ 


ft  provided  with  the  usual  norms  ||  •  ||^  and 
of  those  functions  in  which  vanish  on  the  boundary  in  the  generalized  sense.  The 

inner  product  of  is  denoted  by  (•#•).  Finally,  by  c  we  mean  a  generic  positive 

constant  which  may  vary  with  the  context  but  is  always  independent  of  h. 

Under  assumption  (T) ,  we  have  the  well  known  mean- square-error  estimates 

(1.2) 


lu  -  V^k,2-  ch2"kHul 


2,2 


,  k  =  0,1  , 


and  the  uniform-error  estimates  (see  (4],  [8],  [6],  [1],  17) ) 

(1-3)  Hu  -  V"k,»-  ch2~kln  £  llull2,»'  k  “  0,1  • 

From  (1.2)  and  (1.3)  one  may  conclude  by  an  interpolation  argument  that  for  2  <  p  < 
the  L^  error  behaves  like  (see  (8] ) 


(1.4) 


lu  -  v'lp-  cl,2(ln  £>1_2/PIM 


2,P 
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It  hu  been  considered  as  a  challenge  among  the  specialists  to  remove  the  additional 

logarithm  in  (1.3)  and  (1.4).  This,  in  particular,  since  one  can  show  that  for  higher 

than  second  order  finite  elements  these  estimates  hold  without  the  logarithm  (see  (S), 

®1  2 

(8]).  Also,  for  any  function  u(  n  Wp,  1  <  p  <  «,  the  natural  piecewise  linear 
interpolant  l^u  e  is  well  defined  and  satisfies 


For  the  case  of  linear  finite  elements  Fried  [2]  has  recently  published  an  example  based 
on  radial  symmetry  which  indicates  that  (in  two  and  three  dimensions)  at  .least  the 
pointwise  estimate 

a.®)  ll«  -  v'l«- cl,2ln  jr|lul|2,» 

may  be  of  optimal  order.  However,  this  leaves  the  question  open  whether  the  Lp 
estimate  (1.4)  is  optimal.  In  the  present  paper  we  shall  give  an  answer  to  this  ques¬ 
tion  for  the  model  situation  considered  here  which  is  based  on  the  following  stability 
result: 

•1 

Theorem.  Under  assumption  (T)  the  Rltz  projection  is  stable  in  Wp  for 

2  ^  p  <  »,  namely 

(1-7)  HVII1,P-C||U'I1,P  • 

The  proof  of  the  theorem  will  be  given  in  the  next  two  sections.  One  of  its  consequences 
is  the  following 

•1  2 

Corollary.  Under  assumption  (T) ,  for  any  function  u  *  Wp  ft  Wp  there  holds 

(1.8)  ||u  -  Rhull  l,p  1  chHull2,p'  2<p<-, 

(1.9)  ||u  -  RjU||p  <  cph2IM!2(p-  2  <  p  <  •  . 

Proof,  we  apply  (1.7)  for  u  -  lhu  and  observe  that  a  id  on  Sh  to  obtain 

II  v  xhul'i,p  -  cl'u  ’  *hu|,i,p'  2-p-“  • 

Then,  the  approximation  estimate  (1.5)  implies  (1.8). 

It  prove  (1.9),  we  use  a  duality  argument.  Let  p  *  (2,®),  so  that 
q  m  p/(p-l)  *  (1,2].  On  the  convex  polygonal  domain  0,  the  Laplacian  is  a 


2- 


•1  2 

homeomorphism  from  W  0  w  onto  L  ,  1  <  q  <  2  (see  [3]).  Hence  there  is  a 

q  q  q  — 

°1  2 

v  «  Wq  n  Wq  satisfying 


-Av  "  sgn(u  -  R^u)  |u  -  R^u^  ^  in  ft 


and 

(1.10) 


|v||2#q  lc||Av|jq  -  c||u  -  ^u„p 


P-1 


Using  now  (1.1),  Holder's  inequality, (1. 5) ,  (1.8),  and  (1.10),  we  find 

IP 


V 


(1.11) 


(V(u  -  Rj^u) ,  V(v  -  Ihv)) 
±  Hu-Vlli,pl!v-Ihv>li,q 

1  llu-  Vlll,pChllVl>2,q 


lch2l|u||2,p||v||2tq 
<  ch2  llu | 


IIP"1 


q.e.d. 


'2,p'|u- V'p 

We  remark  on  some  extensions  of  our  results.  The  proof  of  the  theorem  and  to  a 
large  extent  also  that  of  its  corollary  make  use  of  the  fact  that  the  Laplacian  con¬ 
sidered  as  a  mapping 


(1  12) 


“1  2 

w  n  w  +  l 

p  p  p 


is  a  homeomorphism  for  p  <  (1,2  +  a),  where  a  is  some  arbitrarily  small  but 

2 

positive  number.  This  is  certainly  true  on  a  domain  with  smooth  boundary,  say  3ft  «  C  , 
for  all  a  >  0,  and  it  is  known  also  for  convex  polygonal  domains  (see  [31 )  where  a 
depends  on  the  size  of  the  maximum  inner  angle,  u>  <  w.  Our  results  extend  to  more 
general  second-order  elliptic  operators  as  long  as  the  corresponding  mapping  (1.12) 
is  a  homeomorphism.  Ir  the  case  of  a  curved  boundary  the  proofs  become  more  involved 
due  to  the  approximation  of  ft  by  polygonal  domains  ft^.  In  the  case  that  3fi  is 
smooth  one  can  show  that  for  all  p  f  (1,«>]  the  following  refined  estimate  holds: 

,1-1/P  II..  ,1 


(1.3) 


Hv«i.p  -c{||u|li,p,nh  ' h 


)  . 


"  "l.p.JAtib 

From  that  estimate  one  can  again  draw  the  conclusions  (1.8)  and  (1.9),  now  valid  for 
all  p  e  (1,*)  and  p  e  (l,00),  respectively.  The  results  for  1  <  p  <  2  are  proved 


-3- 


via  a  duality  argument  that  makes  use  of  elliptic  regularity  results  that  are  not 
generally  valid  for  non-smooth  boundaries. 

2.  Proof  of  the  Theorem. 

Notation  and  techniques  are  similar  to  those  used  in  [1].  However,  the  key  differ¬ 
ence  is  in  the  type  of  Green's  function  employed.  The  basic  technique  used  in  several 
papers  is  to  reduce  to  the  problem  of  estimating  the  Galerkin  error  g  -  gh  in  approx¬ 
imating  the  solution  of 

-Ag  m  4  in  ft 

where  6  is  the  Dirac  6-function  or  some  approximation  to  it.  The  difficulty  is  that, 
with  piecewise  linear  approximation,  the  error  g  -  gh  contains  a  logarithmic  factor. 

For  example,  it  was  noted  in  18]  that  0  <  c  <_  h  1  (In  h_1) (|g  -  gh  (|  ,  <_  c  as 

1  w1  * 

1 

h  ■*  0.  The  reason  is  that  the  smoothness  of  g  is  such  that  piecewise  linears  fail 
to  afford  optimal  approximation  (whereas  higher  degree  piecewise  polynomials  would 
yield  an  approximation  rate  devoid  of  the  logarithmic  factor) .  The  remedy  here  is  to 
consider  instead  a  "derivative"  Green's  function,  satisfying 

I'  " 

(for  each  i  ■  1,2).  Now  g  is  more  singular,  and  piecewise  linears  afford  optimal 

approximation,  albeit  at  a  slower  rate.  He  now  turn  to  the  details. 

°1 

Let  u  *  Wp,  2  <_  p  <_  «■,  be  given.  We  pick  any  point  z  e  ft  contained  in  the  . 

interior  of  some  triangle  K  «  s.  ,  and  denote  by  8  any  of  the  operators  3/3x. , 

z  n  i 

i  •  1,2.  Because  of  assumption  (T)  there  is  a  function  6  «  C  (K  )  such  that 

z  o  z 

(2.1)  /  6zdx  -  1,  |v*6j  <^ch’2_k,  k  -  0,1,...  , 

where  the  constant  c  does  not  depend  on  z  or  h.  Then,  by  construction, 

(2.2)  *htE)  “  (*h'4»),  V*h*Sh- 
Correspondingly,  we  define  g  e  by 

Z  2 

(2.3)  (Vg2,V*)  -  <«z,3*),  V  *  t  W*  . 


-4- 


Clearly,  gz  is  a  regularized  derivative  of  the  Green's  function  of  the  Laplacian  on 
SI.  Using  this  notation,  we  have 

(2.4)  aR^ufz)  »  07R^u,7gz)  »  (Vu,7R^g^) 

-  (3u,6  )  -  (?u,7(g  -  Rg  ))  . 

z  z  n  z 

We  introduce  the  weight  function 

(2.5)  oz<x)  “  (|x  "  z|2  +  <2h2)1/,2(  f  >  1 

where  the  parameter  <  will  be  chosen  appropriately  large,  «  >  «,>  1,  but  independent 
of  hi  We  note  that  from  now  on  the  generic  constant  c  is  also  independent  of  <  and 
z  e  ft,  and  of  the  parameter  a  e  (0,1]  introduced  below. 

Suppose  temporarily  that  p  <  Applying  Holder's  inequality  to  the  terms  in 
(2.4),  we  obtain  for  any  a  e  (0,1]  that 

Eri 

l<Vu,V(gz  -  Rhgz) )  I  <  (/oz2_a|Vu|Pdx)1/P( /oz2_adx)2p  (  /  <J2+a  |  V (g^  -  Rhgz>  1 2dx) 1/2 

2zl 

<  c  (a'V“) 2p  ^  (  /  o"2*3 1  Vu  |pdx) 1/p  , 

where 

Mb  -  max  (/a2+a|7<gz  -  X^gJ  1 2dx) 1/2  . 

Furthermore , 

Edl 

1  Ou,  6  )  |  <  (/  |3u|pdx)1/p(/  |6  |p_1  dx)  P 

Kz  Kz 

_2 

£  ch  p  (/  |7u|Pdx)1/,p  . 

K 

z 

We  apply  the  above  estimates  with  (2.4),  raise  to  the  p-th  power  and  integrate  with 
respect  to  z  e  ft  to  obtain 


-5- 


||  3R.  u  |1  ^.c(h"2  //  |Vu|Pdx  dz)1/,£>  + 

P  Kz 

P-2 

*  c(aV)2p  Hh(//«r2_a!^lp«*  ^)1/p  . 


Thus,  by  interchanging  integration,  we  find 


(2.6)  ll3RhUllp  -  cH7uilp(1  +  8  1/2h  ' 

where  the  constant  c  is  obviously  independent  of  p.  Estimate  (2.6)  is  also  easily 
seen  to  hold  for  p  »  ®  using  the  above  techniques.  Now,  to  prove  the  assertion  of 
the  theorem,  we  have  to  show  that 

(2.7)  \  -  max  (/  o2+“ |v(gz  -  R^)  |2dx)1/2  <  Caho/  , 

z«  fl 

for  a  proper  choice  of  a  f  (0,1). 

To  prove  (2.7),  we  need  some  preparations.  From  now  on,  we  drop  the  subscript  z 

and  simply  write  c,g  for  °2<gz- 

The  weight  function  a  satisfies 

(2.8)  —  CCrl  11  “  0,1,2,...  . 

Here  denotes  the  tensor  of  *-th  order  derivatives  of  o.  Moreover,  for  k  *i 

sufficiently  large,  one  has  that  (see  [5)) 

(2.9)  max  (max  o  (x)/min  o(x)J 

K*tt.  xtK  X«K 

n 

holds  uniformly  for  z  «  0.  For  any  function  v  e  wj  n  [  n  W2(K)1  the  natural  piece- 

Keir. 

h 

wise  linear  interpolant  Ihv  e  S h  is  well  defined  and  satisfies 

(2.10)  l|V(v-  xhv)||2;K  <  ch||72v||2;K,  n.h. 

Combining  (2.10)  and  (2.9),  one  easily  sees  that  the  following  holds: 


(2.11) 


/  oS|v(v  -  I  v)  |2dx  _<  ch2  /  cB|72v|2dx  , 


where  the  abbreviation  used  is 


/’  ...  dx  -  l  /  ...  dx 


K*Th  K 
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To  prove  (2.7),  we  set  1 1/  «  o2+a(g  -  R^g)  and  we  use  (1.1)  to  obtain 

/  o2+a  | V  (g  -  R^g)  1 2dx  -  /  V  (g  -  R^m*  -  I^Jdx  +  j  f  Ao2+“(g  -  R^dx  . 

Thus, 

(2.12)  /  a2+a  I  V(g  -  R^g)  |  2dx  <_  /  a  2  01 1 V ( —  1^)  |  2dx  +  c  /  aa(g  -  Rhg)2dx  . 

From  (2.11),  we  get  by  a  simple  calculation  that 

/  o  2  a|V(i(i  -  I^tji)  |  2dx  ch2  /  a2+a  |  V2g|  2dx  + 

+  ck  2(f  a2+a|V(g  -  R^g)|2dx  +  /  aa(g  -  Rhg)2dx)  . 

We  insert  this  estimate  into  (2.12)  and  find  that,  for  k  <2  sufficiently  large, 

(2.13)  /  o2+a|v  (g  -  R^g)  |  2dx  <_  ch2  /  o2+a  |  ?2g  |  2dx  +  c  /  aa(g  -  R^g^dx  . 

To  handle  the  second  term  on  the  right  side  of  (2.13),  we  employ  a  duality  argument  in 

°1 

weighted  norms.  For  fixed  h,  let  v  (  be  the  solution  of  the  auxiliary  problem 

(2.14)  -Av  »  oa(g  -  Rj^g)  in  0  . 

2 

Since  (2  is  convex,  it  is  guaranteed  that  v  e  W2>  Moreover,  in  Section  3  we  shall 

°1 

show  that  the  following  weighted  a  priori  estimate  holds  for  all  v  e  (0)  such  that 
Av  f  W2(«) : 

(2.15)  /  o  2  a|V2v|2dx  <_  ca  1(ich)  /  a2  a|VAv|2dx  . 

Consequently, 

(2.16)  /  o'2'“|  V2v|2dx  <_  ca_1(Kh)'2  /  (o2+a|v(g  -  F^g)  | 2  +  0%  -  R^^Jdx  . 

Using  (2.14),  we  have 

/  oa (g  -  Rhg)2dx  =  /  V (g  -  R^gJ-Vtv  -  I^vjdx 

1  (/  o2+a[v(g  -  R^g)  j2d3t)1/,2(/  o  2  a  |  V  (v  -  Ihv)  1 2dx) 1^2  . 

Then,  by  (2.11)  and  (2.16),  choosing  <  ^  <3  yields 

(2.17)  /  0° (g  -  R.g)2dx  <_  c(a<)  1  /  o2+1|v(g  -  R^g)  |2dx  . 


-7- 


We  insert  (2.17)  into  (2.13)  and  choose  again  <  ^  <4(a)  sufficiently  large  to  obtain 

(2.18)  /  c2+<3|V(gz  -  Rhgz)|2dx  _<  ch2  /  c2+Cl  |  V^g  1 2dx  . 

Thus,  we  have  reduced  the  proof  of  (2.7)  to  an  a  priori  estimate  of  the  form 

(2.19)  /  c2+a|'72g j2dx  <_  c^h01  2  . 

This  estimate,  however,  is  an  obvious  consequence  of  the  a  priori  estimate 

(2.20)  /  o2+a|V2g|2dx  £c  /  c2+a(3i)2dx  +  ca  ^(rh)  2  /  o2+a62dx  , 
which  will  be  proven  in  Section  3,  for  0  <  a  ^  sufficiently  small. 

3.  Some  Weighted  ft  Priori  Estimates. 

»1  °1  2  °1 
Let  functions  f  (  and  b  e  [Wj]  be  given,  and  let  v  e  W2  be  such  that 

(3.1)  -Av  =«ft  div  b  in  ft  . 

2  2  1/2 

If  a  =  ( | x  —  z |  +  G  )  '  is  the  weight  function  introduced  in  Section  2,  then  we 

have  the  following 

Lemma.  For  any  convex  polygonal  domain  ft,  there  exists  an  e  (0,1)  such  that 
for  all  parameter  values  a  e  (O.a^)  the  following  a  priori  estimates  hold, 

(i)  if  f  s  0; 

(3.2)  J  o2+a  |  V2v|  2dx  c  J  o2+a|div  b|2dx  +  ca  2  /  o2+a|b|2dx  , 

(ii)  ijE  b  =  0: 

(3.3)  /  a  2  a|v2v|2dx  ^  ca  2  J  o2  a|vf|2dx  . 

Proof.  (i)  To  prove  (3.2),  we  estimate 

J  o  l72V'  dx  —  /  l72  0  V) ]  dx  +  c  J  {o  I 7v I  +o  v  }dx  . 

Since  ft  is  convex,  we  have  the  standard  b^  a  priori  estimate 

IMI2,2  i  ciUwll2'  W  e  w2  n  w2  . 

Applying  this  to  o1+a^2v,  we  find  by  a  simple  calculation  that 
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(3.4)  /  o  a  |  V^v  |  dx  <_  c  j  a  a|div  b|  dx  +  c  /  (o  |  Vv  |  +  au  v  )dx  . 

Furthermore, 

/  aa|Vv|2dx  =  /  VvV(cav)dx  +  j  /  Aoav2dx  , 


and  hence,  using  (3.1), 


(3.5) 

Combining  (3.5)  with  (3.4),  we  arrive  at 


/  oa|Vv|2dx  _<  c  /  o2+ajdiv  b)2dx  +  c  j  0°  2v2dx  . 


(3.6) 


j  a2+a  |  V  v  |  ^dx  <_  c  /  a^+a|div  bj^dx  +  c  / 


2+a  i 


a-2  2 
a  v  dx  . 


Next,  we  apply  Holder's  inequality  to  obtain 


(3.7) 


/  aa  2v2dx  _<  (/  a  A  adx) 


-2-a  \  (2-a)/  (2+a)  n  u2 


(2+a)/a 


-1  -cc  <2-a)/(2+a)  ,,,,2 


<  c  (a  t,  ) 


(2+a)/a 


°1  2 

We  have  already  noted  that  the  Laplacian  is  a  homeomorphism  from  <3  onto 

°1  2 

for  all  q  e  (1,2).  Hence,  there  is  a  we  W(2+ct)/2  n  W(2+a)/2  satlsfyln9 

-Aw  =  sgn  (v)  jv  |  2/>Ci  in  ft  , 


and 

(3.8) 


W  2,(2+a)/2  —  C  (2+a)/2  ‘ 


Then,  we  have  via  Holder's  inequality,  Sobolev's  inequality,  and  (3.8)  that 
(2+a)/a  =  (Vv,Vw)  =  (b,7w)  i  (4+2a)/(2+3a)  HW^1,  (4+2a)/(2-a) 

-  (4+2a)/  (2+3a)  ^W^2,  (2+a)/2  -  (4+2a)/ (2+3a)  (2+a)/a 

Thus,  we  obtain 
(3.9) 


(2+a)/a  -  C'lb'l  (4+2a)/(2+3o) 


Now,  again  by  Holder's  inequality, 

II  bll  (4+2a ) /  (2+3a)  ^  V  o2+“  |  b|  2dx) 1/2  (/  o"  t2+“>  ^x)  a/ (2+“> 

(3.10) 


<  c  C 


-<4+a  )/(4+2 a)g  a2+«  | b|  2^  1/2 


-9- 


Combining  the  estimates  (3.10)- (3. 7) ,  we  obtain  that  for  a  e  (0,1] 


(3.11) 


/■  a- 2  2  -1-2  r 

Jo  v  dx  <_  ca  c  / 


b  dx  . 


This  together  with  (3.6)  proves  the  estimate  (3.2)  for  the  choice  -  1. 

(ii)  To  prove  (3.3),  we  apply  Holder's  inequality  as  follows: 

(3-12)  /  o-2->2v|2d*  -  (/  o-(2+a,/“dx)a||V2v||2/(i-a) 

—  c^a~2  II v  II 2, 2/  (1-a)  . 

*  °1  2 

Above,  we  have  noted  that  the  Laplacian  is  a  homeomorphism  from  w2/(i  a)  n  W2/(l-a) 
onto  h2/ (1  a)  for  a  e  ( 0 ,  ,  where  1  >  >  0  is  determined  by  the  maximum 

inner  angle  of  ft.  Thus,  for  a  e  (O.a^l,  we  have  that 


(3.13) 


Ml2,2/(l-a)  —  c  il  2/  (1-a) 


By  Sobolev's  inequality  combined  with  Poincare's  inequality  (notice  that  Av  e  W2> , 


lAvil2/(l-a>  —  °  ^Av  2/  (2-a)  * 


(3.14) 

We  apply  again  Holder's  inequality  to  obtain 

| ] VAvr  11 2/ ( 2— a)  -  (/  o-(2‘a)/U'a)dx)(1‘“>/2(/  o2'“|7Av|2dx)V2 

(3.15) 

<  ca-1/2C-°/2(/  c2-“|VAv|2dx)1/2  . 

Combining  the  estimates  (3.15)  -  (3.12),  we  finally  reach  the  desired  estimate 

r  -2-a i „  |2,  -1—2  r  2-ai„.  1 2 

)  a  l72v|  dx  <_  ca  C  Jo  |VAv|  dx  , 


valid  for  a  e  (0,a^). 


q.e.d. 
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